Robust quantum gates and a bus architecture for quantum computing with 

rare-earth-ion doped crystals. 
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We present a composite pulse controlled phase gate which together with a bus architecture im- 
proves the feasibility of a recent quantum computing proposal based on rare-earth-ion doped crys- 
tals. Our proposed gate operation is tolerant to variations between ions of coupling strengths, pulse 
lengths, and frequency shifts, and it achieves worst case fidelities above 0.999 with relative variations 
in coupling strength as high as 10% and frequency shifts up to several percent of the resonant Rabi 
frequency of the laser used to implement the gate. We outline an experiment to demonstrate the 
creation and detection of maximally entangled states in the system. 

PACS numbers: 03.67.-a, 33.25.+k, 82.56.Jn 



Introduction 

The Rare Earth Quantum Computer (REQC) recently 
proposed by Kroll and coworkers is a solid state ensemble 
quantum computer based on laser-addressed rare earth 
ions embedded in cryogenically cooled inorganic crystals 
PI @ . Internally, the system consists of a macroscopic 
number of independent instances of the same quantum 
computer. The qubits, represented by ground state hy- 
pcrfine levels in the rare earth ions, are divided into chan- 
nels according to the inhomogeneous frequency shifts of 
an optical transition in the ions. Every instance of the 
quantum computer has one representative from each ac- 
tive channel, and the instances are operated in parallel by 
addressing the channels with resonant optical radiation. 

The inhomogeneous shift of the ions in a given channel 
will not be quite identical, and there will in general be 
fluctuations in their coupling to an applied coherent field. 
Also, the static dipole coupling used to mediate multi 
qubit gate operations will differ between instances. As a 
consequence, we must use gate operations that are stable 
with respect to variations in these parameters which, as 
we will demonstrate in this paper, may be achieved by 
using composite pulses and phase compensating opera- 
tions. 

One of the exciting features of REQC is that the size 
and coupling topology of the quantum computer is not 
defined by the crystal, but rather chosen in an initializa- 
tion stage at each start-up of the system. The choice of 
architecture determines the number of instances available 
in a given crystal, and thus ultimately the scaling prop- 
erties of the system. We propose to choose a bus-based 
architecture, as this will both allow for a higher num- 
ber of instances and simplify the gate implementation as 
compared to the originally proposed architecture. 

We conclude the paper with a proposal for an experi- 
mental demonstration of creation and detection of max- 
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imally entangled states in REQC systems. 



QUANTUM COMPUTING WITH RARE 
EARTH IONS 



Rare earth ions embedded in cryogenic crystals have 
a number of features making them suitable for quantum 
information processing [3|: 

• Ground state hyperfine levels with a lifetime of 
hours and decoherence times up to several ms. We 
will use three such states, labeled |0), |1), and |aux), 
to implement quantum registers and for parking 
unwanted ions. 

• Optical transitions with homogeneous line widths 
on the order of kHz are inhomogeneously broad- 
ened to several GHz, allowing us to address a large 
number of independent channels. 

• The crystal-embedded ions have large static dipole 
moments with interaction energies up to several 
GHz. This interaction is ideal for implementing 
gate operations of the dipole blockade type. 

In the remainder of this section we will briefly intro- 
duce the basic ideas of REQC, as originally described in 
Ref. 0. 



A. Dynamical architecture selection 

The architecture of the REQC system is selected at 
start-up by an initialization procedure. The desired end- 
point of this process is a large number of independent 
instances of the chosen quantum computer, each instance 
being a group of ions with one representative from each 
active channel and couplings between the ions as required 
by the chosen architecture. 
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The initialization proceeds in two steps: channel 
preparation and identification of quantum computer in- 
stances. In both of these steps unwanted ions are deac- 
tivated by transferring them to off- resonant, metastable 
states. 

a. Channel preparation. A channel refers to a large 
number of ions distributed throughout the crystal, all 
having the same inhomogeneous shift and coupling 
strength within the inhomogeneously broadened optical 
transition used to access the ions. The channel prepara- 
tion aims to deactivate all dopant ions close to resonance 
with a given channel and to transfer all members of the 
channel itself to their |0) state. 

This can be achieved by means of spectral hole burn- 
ing techniques, and widths of the final channel structure 
as low as 50 kHz have been obtained experimentally for 
materials similar to those considered for use in REQC 
0. 

b. Instance identification. After a successful initial- 
ization, each ion will only be interacting with ions from 
other channels, allowing us to ignore "excitation hop- 
ping" transitions 0] , as these will not be energy conserv- 
ing. As a consequence we can model the dipole coupling 
as simple couplings between the excited states: 

Vdipole = ^ Yl 9v> (|e) (e|) M ® (|e> (e|)„ , (1) 

where the sum is over all pairs of ions. To be precise 
about the objectives of the instance identification pro- 
cess, we will consider ions fi and v to be coupled if g^ 
exceeds a threshold gt determined by the chosen imple- 
mentation of the gate operation. 

The goal of the instance identification procedure is to 
transfer ions, which are in an active channel but not 
members of a valid instance, to their auxiliary state |aux). 
One way to achieve this is to go through the following 
procedure for each pair, of channels required to be 

coupled: 

By applying a 7r-pulse to ions in channel i we trans- 
fer the |0) population to the |e) state, thus shifting the 
excited state energy of all ions coupled to a channel i 
ion. By means of a frequency sweep or a comb of n ro- 
tations, all channel j ions which are shifted less than gt 
are now transfered to their excited state |e), after which 
the channel i ions are returned to |0). We now wait for 
the excited channel j ions to decay, which will transfer 
part of the ions to the inactive |aux) state. 

By repeated application of this pulse sequence, we can 
deactivate an arbitrarily high fraction of the channel j 
ions which are not coupled to a channel i ion. After this 
has been achieved, we repeat the process with the roles 
of channel i and j interchanged, and afterwards proceed 
to apply the same procedure to all other edges of the 
coupling graph to finally arrive at the desired initialized 
REQC system. 



B. Gate operation 

In general, the coupling strengths, g^ Vl will differ be- 
tween instances, requiring us to use gate operations that 
do not depend on the precise magnitude of the coupling 
strength. 

One gate operation with this quality is the controlled 
phase shift based on the dipole blockade effect |5|. As- 
suming all ions not participating in the operation to be in 
their qubit states, |0) and |1), and thus decoupled from 
the operation, we can implement a controlled phase shift 
in its simplest form by the following pulse sequence: 

P$(n,K)p£>(0,2n)P$(0,K), (2) 

with pf$ (4>, 6) representing the effect of a resonant pulse 
of area and phase cj> applied on the \a) — \b) transition 
of ions in channel i. 

For two coupled ions /i and v, residing in channels i and 
j respectively, the effect of performing the pulse sequence 
(|2J) would be the following: If ion /i is initially in the 

0}^ state, the P$ (0, 7r) pulse transfers the ion to the 

excited state |e) and thus shifts the P^(0, 2ir) pulse out 
of resonance, causing the system to return to the initial 
state after the last 7T pulse. If, on the other hand, ion fi 
is initially in the |1) state, it is not transferred to |e) 
and the 2ir pulse is resonant and causes a tt phase shift 
on the state. The effect of the full gate operation 
on the qubit space is consequently a n phase shift on the 
11) state: U C ps = 1 — 2 [11> (11|. 

II. HIGH FIDELITY GATE OPERATIONS 

Gate operations for the REQC system face a number 
of challenges due to the fact that they operate simulta- 
neously on a number of not quite identical instances of 
a quantum computer: Due to the finite channel width, 
ion pL will in general be detuned by a small amount 5^' 
from the central channel frequency. Furthermore, the ex- 
perienced Rabi frequency, fl^, will differ slightly from 
the average Rabi frequency f2o due to laser field inhomo- 
geneities and local variations in dipole moments. 

In this section we will show that by taking advantage 
of the fact that and Qq^ /VLq are constant in time 
for each ion, we can design pulse sequences that perform 
almost the same operation on each instance. 

A. Composite rotations 

The pulse P ie (</>, 6) is driven by a Hamiltonian Hi = 
i$l • er( ie ) with cr( ie ) signifying the Pauli-matrices in the 
{\i) , |e)} basis and fl = Qofi^, where is a unit vector 
in the x — y plane with azimuthal angle <f>. 

To apply the pulse P ie {4>, 0) we engage the field for a 
period so that an ideal reference ion, £, with <Jw = 
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and ft, 



(?) 



Qo will be rotated by an angle 8 around 



fifj, as desired. In general, however, the ions will react 
differently to the pulse due to their different detunings 
and coupling strengths. 

The problem of taking all the ions through the same 
evolution when they react differently to the pulses has 
been studied in great detail in the magnetic resonance 
community [||. Inspired by the discussion in Ref. we 
have used the BB1 pulse sequence to replace a single 
pulse -P(0, 9) with the following sequence of pulses: 

IWM) = 

P(0, 6/2) Pfc, tt) P(30 c , 2tt) P(<t> c , tt) P(0, 0/2). (3) 

For our reference ion, £, the unitary evolution ^Pbbi (</>, 0) 
caused by the Pbbi(0, composite pulse is seen to be 
exactly identical to the evolution caused by P(<fi, 0). The 
use of five pulses for this simple task is justified, however, 
if we instead consider the evolution ^Pbbi of a general 
ion subject to the Hamiltonian 
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(m) 



-<5 (ai) |e) (e 



5°° 



n</> 



r (ie) 



(4) 



In this case we find that with the optimal value, <p c — 
± cos -1 (— 9/4n), m Pbbi(0, 0) is almost constant over 
a large range of values of $W /Q and fiW/fl , while 
M P(</>, 6*) changes quite rapidly. 



B. Robust gate operation 



To asses the gate performance we have compared the ef- 
fect of the gate to the desired gate operation, Uqps, in 
terms of the worst case fidelity, defined as the minimal 
overlap between the actual outcome of the pulse sequence 
and the desired outcome of the gate operation: 



•^(t/cPs/P) = min 
Wen 



U, 



CPS 



•P 



(6) 



Since we know that the starting point of the gate op- 
eration will be a superposition of the ground hyperfine 
states, we have not minimized the expression © over the 
full Hilbert space, but rather restricted TC to the qubit 
space. Note that this modification ensures that any pop- 
ulation in the excited state after the gate operation is 
counted as a loss of fidelity as it should be. The com- 
putation of the fidelity is discussed in more detail in the 
appendix. 

As we see from Fig.^the pulse sequence Pops obtains 
high fidelities over a much larger parameter space than 
the simple gate operation described by Eq. This is 
highly desirable, as the minimal fidelity among the in- 
cluded instances determine the scale-up needed to per- 
form error correction Q. Not too surprisingly, the sen- 
sitivity to variations in fl is improved the most, as this 
is the type of error best dealt with by the BB1 pulse se- 
quence. For realistic parameters of the REQC system, 
a reduced sensitivity to 5 variation would be more use- 
ful; whether this can be achieved by means of composite 
pulses is a point of further study. 



For the two-level Rabi problem there is a global phase 
factor depending on the detuning which plays no observ- 
able role. In our three-level system, however, this phase 
will lead to a dephasing between the qubit level \i) cou- 
pled to e) and the other qubit level. To compensate 
this, we must symmetrize the desired pulse sequence in 
a suitable way, to allow both levels to pick up the same, 
unknown, phase contributions. 

In the case of the controlled phase shift @, we have 
arrived at the following symmetrized version: 

P« (tt, tt) P« (0, tt) Pg (tt, tt) Pg (0, tt) 

P«(0,7t)P <:W) 
P « (tt, tt) P$ (0, tt) P«> (0, tt) P« (0, n) 

P$(0,n). (5) 

(i i) 

For the reference ion £ the P^pg pulse sequence is seen 

to be equivalent to P ( * } (tt, tt) P^ } (0, 2tt) P$ (0, tt) , which 
is exactly the basic controlled phase shift operation J2J), 
but we expect it to perform better for a general ion. 

(i i) 

Implementing all the pulses of P^pg by composite BB1 
pulses , we do indeed obtain a very robust implemen- 
tation of the controlled phase shift as illustrated in Fig.^ 



III. THE BUS ARCHITECTURE 

As the architecture of an REQC system can be chosen 
at will, the question remains of which architecture to 
choose. 

The fully interconnected "cluster" architecture sug- 
gested in the original REQC proposal of course has the 
minimal topological distance between qubits. On the 
other hand, a star topology with one central qubit cou- 
pled to the n — 1 remaining qubits, as illustrated in 
Fig. |21 would reduce the number of required couplings 
from n(n — l)/2 to n — 1, thus increasing the number 
of available instances in a given crystal, while still main- 
taining a topological distance of only 2. 

Since the outer qubits in the star topology are not di- 
rectly coupled, two-qubit gates between those must be 
mediated by the central qubit acting as a bus. To be 
specific, a bus-mediated controlled not gate can be con- 
structed as 



Bus 




with the vertical lines signifying controlled phase shifts. 

In addition to the better scaling properties of the bus 
architecture, its main advantage is that the bus qubit is 
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FIG. 1: Calculated worst case fidelities @ of two implementations of the controlled phase shift: a) The simple implementation 
J2J), and b) the Pops pulse sequence @. The fidelity is plotted as a function of 8^ = <5' 2 ' and Q : 



(i) 



n 



(2) 



both relative to 



fio, and with g±2 = 100 f2o- Note the difference between the fi-axis limits of the two plots. It is clear from the plots that Pops 
achieves a high fidelity over a much larger parameter space. In particular, Pops is much less sensitive to variations in SI, while 
the sensitivity to variations in 8 does not seem to be significantly improved. 




FIG. 2: Two possible coupling topologies for REQC systems: 
a) cluster topology and b) star topology 



a participant of all multi-qubit gates. This fact can be 
used to ease or improve the implementation of such gates: 
As an example, four times as many pulses are needed on 
channel j as on channel i with the proposed implemen- 
tation of the controlled phase shift (JSJ. If channel j is 
chosen as the bus channel, a dedicated laser system can 
speed up the application of these pulses as compared to 
a tunable laser system able to address any channel. 



IV. PREPARATION AND DETECTION OF 
MAXIMALLY ENTANGLED STATES 

To demonstrate the viability of the REQC concept, 
and in particular the bus architecture, we propose to 
perform an experimental preparation and detection of a 
maximally entangled state. 

We will use an REQC system with star topology: one 
central qubit coupled to n — 1 outer qubits. Starting 
with all n qubits in their |0) state, we apply a composite 
pulse Hadamard operation to the central qubit followed 
by controlled not operations on all the outer qubits con- 
trolled by the central qubit, thus transferring the system 



to the maximally entangled state 

l*Q> = 73(|0> n + |l> w ), ( g ) 

which corresponds to a superposition of the total pseu- 
dospin pointing straight up and straight down. 

The following algorithm for detecting a population of 
the cat state j^o) is very similar to the method used 
by the group of D. Wineland to detect a maximally en- 
tangled state of four ions in a linear Paul trap [tj: By 
rotating the state |^o) through an angle (j> around the 
2-axis we accumulate different phases on the pseudospin 
components: = ^= (|0)" + e"^ n |1)™). An addi- 

tional rotation by 7r/2 around the y-axis now yields a 
state 1^2) with an expected parity, P — TLi(a z )i, given 
by 

(* 2 |P|*2) =cos(n0), 

the detection of the n<j) dependency thus signifying that 
the maximally entangled state has been populated jlCj . 

In a single-instance quantum computing system, such 
as the ion trap setup used in Ref. |9|, we could measure 
the expectation value of the parity as a statistical average 
over many repetitions of the procedure described above: 
after each run we could simply measure the state of each 
qubit, and subsequently compute the parity. Since mea- 
surements in the REQC system yields an ensemble av- 
erage, this approach would not be applicable here: we 
cannot find the expectation value of the parity from the 
ensemble averages of the single qubit parities, ((cr z )j), 
which are as inspection shows. 

Instead we let the bus qubit acquire the parity unitar- 
ily: by sequentially applying controlled not operations 
from each outer qubit to the central qubit we make the 
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central qubit end up in the |1) state in the case of odd 
parity and in the |0) state in the case of even parity. Af- 
ter this, the ensemble average of the bus qubit population 
yields the expectation value of the parity. 

As this section illustrates, readout from an ensemble 
quantum computer is conceptually somewhat more com- 
plicated than readout from a single quantum computer. 
It is worth noting, however, that unlike many other en- 
semble quantum computing proposals, REQC instances 
all start in the same pure state: if we successfully em- 
ploy error correction during a computation all instances 
will end up in the same pure state, allowing us to read 
out the ensemble averages with high signal to noise ratio. 
Perhaps surprisingly, the readout can almost always be 
performed by tricks similar to those employed to detect 
the maximally entangled state: Ensemble quantum com- 
puting is almost as powerful as general quantum com- 
puting. In particular, all problems which may be ex- 
pressed in terms of the hidden subgroup problem (such 
as Shor's factoring algorithm) can be solved using an en 
semble quantum computer 



n) c 

m 



channel populations in their |aux) state), the selection 
could theoretically be performed unitarily. 



APPENDIX: FIDELITY OF UNITARY 
OPERATIONS 

We wish to compare unitary operators U and Uq, by 
determining how closely UqU resembles the identity on 
the Hilbert space Tt. This can be expressed in terms of 
the worst case fidelity: 



T(U ,U) = mux 



ulu 



1> 



(A.l) 



The fidelity can be computed as follows: UqU is uni- 
tary and can consequently be formally diagonalized with 
eigenvalues e l ^ j , j = 1, . . . , n so that < <j>i < ... < 
4>n < 27r. Introducing the maximal eigenvalue phase dis- 
tance A<j> max = max({0j-0j_i}j =2l ..., 7l U{27r+0i-0n}) ) 
the fidelity over TL is given as 



V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have shown that, in the absence of 
decay and decoherence, it is possible to implement robust 
high-fidelity gates for the REQC system. Specifically, the 
phase compensated controlled phase gate based on com- 
posite pulses ©, achieves worst case gate fidelities above 
0.999, even with the coupling strength varying up to 10% 
between instances and channel widths of several percent 
of the Rabi frequency of the field used to manipulate the 
system. Furthermore, we have pointed out that using a 
bus based architecture will simplify implementation by 
allowing the use of an asymmetric laser setup. 

The number of instances of a bus based REQC system 
scales as p n where n is the number of qubits per instance 
and p is the probability of a random ion being coupled 
to a member of a given channel. In the regime currently 
being investigated experimentally, p is several orders of 
magnitude less than 1 . The value of p is affected by gt and 
channel width, which is why we have to use robust gates 
rather than narrow channels and high threshold coupling 
strengths. Higher values of p could be obtained by in- 
creasing the ion density, which would, however, cause a 
decrease in coherence times. By using structured dop- 
ing techniques it might be possible to obtain a higher 
effective p without this adverse effect. Another approach 
to obtaining higher effective p would be to use multiple 
channels for each qubit by guaranteeing each instance to 
have exactly one member ion from a group of channels 
assigned to each qubit. 

The instance identification protocol described in Sec.[l] 
could be made much more efficient: Since the system 
starts in a pure state (all ions in the channels in their |0) 
state), and also ends in a pure state (all instance mem- 
bers in their |0) state, and all other ions from the initial 



HUo,u) 



cos 2 (A0 max /2) if A</> max > tt, 







otherwise. 



(A.2) 



To see this, we expand the state vector \ip) on the 
eigenbasis {|j}} of U%U: \ip) = J2j c j Thc fidelity 
then takes the form 



T(U ,U) =min 



(A.3) 



with the minimum taken over all non- negative pj — \cj\ , 
so that V pj = 1. 

Eq. (|A.3(1 allows us to interpret the fidelity geo- 
metrically in the complex plane, as the set of points 
Pj ex P(^j)} form a convex polygon with vertices 
in the eigenvalues {e l * J } on the unit circle. The fi- 
delity corresponds to the square of the minimal distance 
from to this polygon. If the polygon is constrained 



to one half-plane, this will be |e 



itp 



a i(4>+A<j> n 



74 



cos (A</> max /2). If the polygon is not restricted to one 
half-plane, it will cover the origin, and the fidelity will 
be 0. 

Note that this method relies on the minimization being 
performed on the whole Hilbert space. If this is not the 
case the method is not applicable, and in Sec. Ill Bl where 
the minimization is carried out over a subspace of the full 
Hilbert space, we have resorted to a numerical search. 
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